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Abstract. A comprehensive temporal analysis has been performed on the 319 brightest GRBs with Tgo>2s 
from the BATSE current catalog. The GRBs were denoised using wavelets and subjected to an automatic pulse 
selection algorithm as an objective way of identifying pulses and quantifying the effects of neighbouring pulses. 

The number of statistically significant pulses selected from the sample was greater than 3000. The rise times, 
fall times, full-widths at half maximum (FWHM), pulse amplitudes and pulse areas were measured and the 
frequency distributions are presented here. All are consistent with lognormal distributions provided the pulses are 
well separated. The distribution of time intervals between pulses is not random but compatible with a lognormal 
distribution when allowance was made for the 64 ms time resolution and a small excess (5%) of long duration 
intervals that is often referred to as a Pareto-Levy tail. The time intervals between pulses are most important 
because they may be an almost direct measure of the activity in the central engine. Lognormal distributions of 
time intervals also occur in pulsars and SGR sources and therefore provide indirect evidence that the time intervals 
between pulses in GRBs are also generated by rotation powered systems with super-strong magnetic fields. 

A range of correlations are presented on pulse and burst properties. The rise and fall times, FWHM and area 
of the pulses are highly correlated with each other. The pulse amplitudes are anticorrelated with the FWHM. 

The time intervals between pulses and pulse amplitudes of neighbouring pulses are correlated with each other. 

It was also found that the number of pulses, N, in GRBs is strongly correlated with the fluence and duration 
and that can explain the well known correlation between duration and fluence. The GRBs were sorted into three 
categories based on N i.e. 3<N<12, 13<N<24 and N>25. The properties of pulses before and after the strongest 
pulse were compared for three categories of bursts. No major differences were found between the distributions of 
the pulse properties before and after the strongest pulse in the GRB. However there is a strong trend for pulses to 
have slower rise times and faster fall times in the first half of the burst and this pattern is strongest for category 
N. This analysis revealed that the GRBs with large numbers of pulses have narrower and faster pulses and also 
larger fluences, longer durations and higher hardness ratios than the GRBs with smaller numbers of pulses. These 
results may be explained by either homogeneous or inhomogeneous jet models of GRBs. The GRBs with larger 
number of pulses are closer to the axis if F varies with the opening angle of the jet and the imprint of the jet is 
preserved in the pulse structure of the burst. The distribution of the number of pulses per GRB broadly reflects 
the beaming by the jet. 
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1. Introduction 


Much of the recent progress in the study of gamma-ray 
bursts (GRBs) results from the detection of bursts with 
good location accuracy by BeppoSAX that enabled the 
detection of counterparts at other wavelengths. The sub¬ 
sequent redshift determination of bursts have established 
that these bursts are at cosmological distances (posta 


of the star formation history of the universe (Lamb & 
Reiclrart, 2000: Hanlon et al., 2000; Berger et al., 2001). 


However not many redshifts are known and there is still 
much work to be done to determine the mechanisms that 
produce these enigmatic events. 

The most plausible GRB progenitors are expected to 
be a newly formed black hole (BH) surrounded by a 


et al., 1997 ; van Paradijs et al., 1997 ). GRBs seem to be temporary accretion disk ( Rees, 1999 ; Meszaros, 2001 


connected to massive stars and become powerful probes pastro-Tirado, 200l[) . The most popular models include 
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the merger of a neutron star (NS) and a NS (Eichler et al. 


1989 ; R.uffert fe Janka, 1999| ), NS and a BH ( Paczynski 
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1991), BH white dwarf merger ( foyer et al., 1999|) and 


models of failed supernovae or collapsars (MacFadyen & 
Woosley, 1999; Paczynski, 1998|). An important exception 


is the model in which the GRB energy is provided by a 


newly formed neutron star (Usov, 1992; Thompson, 1994). 
Various explanations have been put forward for the com¬ 
plicated structure of the light curves. These range from 
internal shocks, caused by variations in the velocity of the 
outflow ( |Rees fe Meszaros, 1994 ; Piran, 1999| ), to external 
shocks, caused by interactions with an external medium 


found an anticorrelation between the time delay in the 
arrival times of hard and soft photons in pulses and the 
luminosity of the GRB. 

The light curves of GRBs are irregular and complex. 
Statistical studies are necessary to characterise their prop¬ 
erties and hence to identify the physical properties of 
the emission mechanism. The statistical methods used for 
temporal studies can be broadly divided into four cat¬ 
egories: (1) fits to individual puls es in the GRB using a 
number of pu lse shape parameters ( Norris et ah, 1996 ; Lee 


( Mcsziros fe Rees, 1993 ; Dermer fe Mitman, 1999). In the et al., 2000fc| Ja|) ; (2) a non-parametric approach to pulse 


internal shock model the instabilities in the wind leads to 
shocks which convert a fraction of the bulk kinetic energy 
to internal energy remote from the central engine. A tur¬ 
bulent magnetic held then accelerates electrons which ra¬ 
diate by synchrotron emission and inverse Compton scat¬ 
tering, generating the GRB. Many of the observed features 


shapes in GRBs ([McBreen et al., 1994 ; Hurley et al., 1998 
Young et al., 1995| ; fouilligan et al., 1999 ); (3) the average 
s tatistical properties of G RBs using a peak-aligned profile 
( Stern fe Svensson, 1996|); and (4) the average power spec¬ 
tral density of GRBs ( Belli, 1992 ; |Beloborodov et al., 2000 


in bursts can be reproduced in the internal shock models of et al., 1994) revealed that lognormal distributions can ade- 


:00l|). 


|Chang fe Yi, 2000|) . One of the first studies ( McBreen 


GRBs (fori & Piran, 1997: 

Kobayashi et al., 1997 

Daigne 

quately describe the properties of GRBs. Subsequent stud- 

& Mo< 

hkovitch, 1998; 

Panaitescu et al., 199S; 

Downes 

ies ( 

Li & Fenimore, 1996; 

Hurley et al., 1998 

; iQuilligan 


et al., 

A variety of analytical techniques has been applied to 
the temporal and spectral profiles of GRBs which place 
constraints on the observed distributions which models 
must satisfy. The impressive results from these studies in¬ 
clude (1) hard to soft evolution ( |Golenetskii et al., 1983 ; 
Borgonovo & Ryde, 2001); (2) the duration-hardness an- 
ticorrelation ( Kouveliotou et al., 1993|); (3) the ternpo- 
ral asymmetry of pulses in GRBs (Nemiroff et al., 1993; 
Link & Epstein, 1996); (4) a bimodal duration distri¬ 
bution of GRBs consistent with two lognormal distribu¬ 


et al., 1999) have confirmed the applicability of lognor¬ 


mal distributions in accounting for the wide range in the 
observed properties of pulses in GRBs. This result is not 
surprising because lognormal distributions arise from the 
product of probabilities of a combination of independent 
events and such conditions apply to the pulse generation 
process in GRBs. _ 


tions (Kouveliotou et al., 1993; McBreen et al., 1994); 


(5) th j_discovgiy__o^_tw3__different types of pulses in 
GRBs( Pendleton et al., 1997|); (6) a cor relation between 
Epeak and intensity) Mallozzi et al., 1995|); (7) energ y de¬ 
pendence of the pulse duration ( Norris et al., 1996 ); (8) 
a relationship betw een the pulse peak energy, E pea k, and 


In a different approach (Beloborodov et al., 2000) used 
Fourier analysis to study the power spectral density of long 
GRBs. This approach revealed that the diversity of GRBs 
is due to random realisations of the same process which is 
self-similar over a range of time scales ([Stern &; Svensson 


1996). The slope of the PSD was -5/3 suggesting that 


the photon fluence ( Liang fe Kargatis, 1996 ; Crider et al. 
19991)71(9) lognormal pulse shapes and tim e intervals be- 


tween pulses in long ( McBreen et al., 1994 ; Hurley et al. 
1998()lnd short GRBs ( McBreen et al., 200l|); ( 10) spectra 


well fit with a Band function ( Band et al., 1993); (11) spec- 
tral hardening before a count rate increase (Bhat et al., 
1994|)~ (12) an x-ray excess in GRB spectra ( Strohmaycr 
et al., 1998|); (13) a correlat ion between complexity and 
brightness ( |5tern et al., 1999 ) and (14) the unique proper¬ 
ties of the pulses and power law relationships between the 


pulse properties and durations of GRBs (McBreen et al. 

2002 ( 7 ] 

While GRBs display hard to soft spectral evolution, 
there is remarkable constancy of the pulses in GRBs 


throughout the burst (Ramirez-Ruiz & Fenimore, 200C; 
Quilligan et al., 1999| ). The temporal and spectral proper¬ 
ties of a few GRBs with known redshift have yielded two 
important results to suggest that GRB properties may be 
related to their luminosities. Ramiriz-Ruiz and Fenimore 
(1999) have shown that more r apidly variable burst s have 
higher absolute luminosities. Norris et al. (2000) have 


GRBs are related to fully developed turbulence. The two 
different approaches are quite similar because the lognor¬ 
mal approa ch has been used to describe fully developed 
turbulence ( Arneodo et ah, 1999 ). 

T he work presented he re expands on the earlier anal¬ 
ysis ( Quilligan et al., 2000| ) and provides new insight into 
the mechanism which generates GRBs. The aim is to pro¬ 
vide a comprehensive description and understanding of 
the pulse properties in GRBs and combine it with other 
studies of the spectral properties. The wavelet analysis 
and the pulse selection algorithm are described in Sect. 
2. The method for comparing the properties of the pulses 
before and after the strongest pulse in the GRB is also 
described in Sect. 2. The results are presented in Sect. 3, 
and discussed in Sect. 4. The conclusions are presented in 
Sect. 5. 

2. Data Preparation 

The dataset used was taken from the BATSE cur¬ 
rent catalogue. The ‘discsc’ files are available a t 


bttp://www.batse.msfc.nasa.gov/batse/grb/catalog/4b/ 


( Paciesas et al., 1999 ). The files contain the data from 
the four energy channels, which were combined into a 
single channel to maximise the signal to noise ratio. The 
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Fig. 1 . Illustration of the background subtraction algo¬ 
rithm. The shaded region indicates the bursting phase of 
the GRB. 


in the original signal. This iteration with scaling of the 
convolved function allows the identification of structure 
with a variety of extents both in the spatial and frequency 
domains. 

If the wavelet function is written as il>(x) then, more 
formally, the transform at a particular scale s, can be writ¬ 
ten as 

W 8 f(s,x) = f*ip a (x) 


shapes of GRB pulses vary little with energy and pulses 
in different energy channels can be added together and 
nearly retain their initial shape. A subset of the BATSE 
catalo gue w as selected based on the criteria ( Norris 
et ah, 1996| ) that the GRB duration was greater than 
two seconds ( Tg 0 > 2 s) and the peak flux P 256 ms > 3.28 
photons cm _2 s _1 . In this way a sample of 324 bursts 
with good signal to noise and clearly resolved features 
was obtained. Five of these bursts could not be analysed 
properly due to data gaps, and so our final sample 
consisted of 319 GRBs. All 319 GRBs were used for 
the timing analyses. The 250 GRBs that were summed 
over only two LAD detectors were used for all analyses 
involving pulse amplitude and area. 


2.1. Background subtraction 

The first step in the data preparation involved select¬ 
ing the appropriate background for subtraction from the 
GRB. The start and end times for each burst were iden¬ 
tified. A further margin of 10 seconds was added to both 
the beginning and end of this chosen section. Two back¬ 
ground sections of duration 30 s were then selected, one 
finishing 20 s before the start of the section containing the 
burst and the other starting 20 s after the end of the burst 
(Fig. 1). These two regions were used to fit a linear back¬ 
ground that was subtracted from the burst section. 


2.2. Denoising technique 

One of the difficulties in analysing the time profiles of 
GRBs is in overcoming the limitations imposed by the 


where ip s (x) = (1 /s)i/j(x/s) represents the wavelet dilated 
by a scale factor s. At each scale 2?, a discrete wavelet 
transform that we denote by W.fj can be computed. For 
the particular choice of scale, s = 2 °, j = 1... J the se¬ 
quence 


{si,, ("$/),<*,} 


presence of noise in the signal and the overlap of the in¬ 


dividual pulses. The transient nature of GRBs also means 
that the usual assumptions for Fourier transform tech¬ 
niques do not hold (Chang & Yi, 2000; |5uzuki et al. 


is called the discrete dyadic wavelet transform of the input 
discrete signal D = (Sif(n)) n£ z- The W^jf components 
provide the detail at each scale, meaning the response of 
the (scaled) wavelet function to the detailed variation of 
the signal. The coarse signal, Sljj. provides the low fre¬ 
quency (slowly varying) component of the signal remain¬ 
ing at scales larger than 2 J . The higher frequency compo¬ 
nents can all be recovered from the dyadic wavelet trans¬ 
form (W* /) 1<7< j between scales 2 1 and 2 J . A fast algo¬ 
rithm for calculating the wavelet transform of a signal was 
developed by |Mallat fc Zhong (1992| ) and implementation 
of this algorithm is at the core of the denoising procedure 
used. In Fig. 2, GRB 920513 is shown together with its 
dyadic wavelet transform for scales s = 2 2 " 4 (the low¬ 
est scale, s = 2 1 , is not shown because it is dominated 
by noise). The transforms at each of the scales are shown 
(Fig. 2c, cl and e), along with the low frequency signal 
S ',24 containing the remaining information for scales i > 4 
(Fig. 2f). 

The function ip(x) used in the algorithm was chosen 
so that the wavelet acts like an edge detector with a delta 
response function to a step edge. In fact Wiy / is propor¬ 
tional to the derivative of the original signal smoothed at 
the scale 2 J . Thus calculating the positions of the modulus 
maxima of the transform, \W$jf\, is analogous to locat¬ 
ing the sharp variations in the original signal. Mallat & 
Zhong (1992|) also developed an algorithm for allowing the 


|200l|). An alternative method of filtering the signal is with 


wavelets. Wavelet analysis was pioneered by Daubechies 
|(1992| ) and others during the 1980’s ( Meyer, 1993| ). 


Wavelets are specific functions that, when convolved 
with the signal under investigation, produce a transformed 
signal that represents the location and strength of varia¬ 
tions within the original data. The convolution is applied 
repeatedly to the data over a range of scales with the 
convolution function gradually stretched to coarser and 
coarser scales, revealing variations at corresponding scales 


reconstruction of a signal given just the modulus maxima 
of the wavelet transforms across a set of scales along with 
the low frequency signal remaining at the coarsest scale. 
This reduced representation can reconstruct an accurate 
copy of the original, using an iterative algorithm that con¬ 
verges quite quickly to acceptable levels. 

A technique for identifying white noise and remov¬ 
ing it without losing any other information was intro¬ 
duced by Mallat fe Hwang (1992| ). This technique relies 
on characterising the behaviour of noise across the var¬ 
ious scales in the wavelet transform using just the in¬ 
formation present in the extrema wavelet representation 
described in the previous paragraph. The change in the 
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amplitude of the extrema between scales allows the num¬ 
ber called the Lipschitz exponent to be calculated. More 
precisely, each extremum describes a particular curve in 
(log(s), log \Wf(s, x)|) space representing its increase or 
decay on all the scales for which Wf(s,x) has been com¬ 
puted. Then the Lipschitz exponent, a, is just the maxi¬ 
mum slope of a straight line that remains above this curve 
(Young et al., 1995). Using results from an analysis of 
white noise (a < — 5 , Mallat fe Hwang (1992| )), as well as 
the investigation described below, the characteristic dis¬ 
tribution of Lipschitz exponents for the noise present in 
BATSE GRB profiles was determined. In general, noise is 
expected to have a negative Lipschitz exponent indicating 
a decreasing amplitude with increasing scale. Extrema in 
the wavelet transform which fall in the range expected for 
noise can be removed using an algorithm developed for 
this purpose. 


As an example consider the extremum at around 110 s 
on the top (finest) scale of GRB 920513 (Fig. 2c). On 
the next scale, the amplitude of this extremum is con¬ 
siderably smaller and by the following scale it is hardly 
detectable. It is clear that this extremum corresponds to 
a noise feature in the original signal and visual examina¬ 
tion of the transform indicates that the amplitude of the 
wavelet transform decays quickly when moving to progres¬ 
sively coarser scales, as expected for noise. The extremum 
at around 95 s corresponds to a pulse with intensity over 
25,000 counts per 64 ms. The amplitude of the wavelet 
maximum increases on coarser scales, contrary to the be¬ 
haviour expected for noise. 


The extrema removal algorithm was combined with a 
simple thresholding procedure based on the analysis of a 
section of the burst background. For each burst being pre¬ 
pared, a section of the same burst well separated from the 
signal was used to generate the wavelet transform of the 
background. The standard deviation of the transform at 
the scale 2 2 (the finest scale not dominated by noise) was 
combined with a significance level (a t) to calculate an am¬ 
plitude threshold for the wavelet coefficients of the signal 
at the same scale. Any extrema in the signal with ampli¬ 
tudes less than this threshold on scale 2 2 were assumed to 
be due to noise and were removed. 


A denoised signal was then re constructed (e.g. Fig. 2a ) 
using the algorithm described by Mallat fe Zhong (1992| ) . 
The reconstructed functions have no spurious oscillations 
or sharp variations, and are very similar to the original 
signal. The residuals between the background subtracted 
burst and the reconstructed burst are shown in Fig. 2b. 
The residuals were calculated assuming that the errors 
in the raw data were Poisson distributed. The variation 
of the residuals in the section containing signal do not 
differ significantly from a section due to background noise 
(Fig. 2b). 


2.3. Pulse selection 

Each pulse was examined to find minima on each side 
which were separated in amplitude from the maximum 
by more than a chosen significance level. If the search for 
minima failed on a pulse it was rejected and the search for 
pulses continued. The algorithm was designed so that the 
maximum and minima finally selected for a given pulse 
were the extreme values in that region of the signal. The 
significance level of the maxima and minima was calcu¬ 
lated by multiplying the error on the counts by a scale fac¬ 
tor, T a . Given two data points n\ and n 2 , where n\ > 712 , 
the points overlap if 


(ni - T CT Ani) < (n 2 + r a An 2 ) 


The errors on the counts in the cleaned signal were 
assumed to be Poissonian. 

The pulse selection process is illustrated in Fig. 3. 
Starting with pulse 2, minima 1 and 3 were easily iden¬ 
tified. Next, pulse 4 was considered. Pulse 4 overlaps the 
adjacent minimum 3 (based on appropriately sized error 
bars) and both turning points were rejected and minimum 
5 was considered as the true minimum associated with 
pulse 2. A minimum, maximum, minimum triplet consist¬ 
ing of points 1, 2 and 5 was found. At the conclusion of the 
analysis the overlapping maxima/minima pairs consisting 
of points 6 and 7, 9 and 10 and 13 and 14 were rejected. 
The algorithm selected points 2, 8 and 12 as the maxima 
and identified associated minima 1, 5, 11 and 15. 

The analysis was then extended to allow the iden¬ 
tification of pulses that were well separated from their 
neighbours such that the overlap from the surrounding 
signal did not significantly affect the profile of the pulse. 
The fraction of the total height, H b , (from pulse to back¬ 
ground) which was above the higher minimum was deter¬ 
mined and a threshold was applied to this fraction, above 
which pulses were considered isolated. In Fig. 3 the pulses 
8 and 12 are not very well separated from their neighbours 
whereas pulse 2 is effectively isolated and not strongly 
influenced by surrounding signal. The fractional isolated 
height of pulse 8 (to which the threshold is applied) was 
obtained using (Fig. 3): 

fi = min (H u H r )/H b = H r /H b = 8/21 ~ 0.38, 
where Hi is the height difference between the pulse 8 and 
the point 5 on the left and H r is the height difference 
between the pulse 8 and the point 11 on the right. 

For pulse 2 (with point 1 being the higher minimum 
the smaller height difference is on the left) the estimate 
is: 

fi = min (H h H r ) = H t /H b = 15/17 » 0.88. 

Thus, if the threshold, 77, were set above ~ 40% then 
pulse 2 would be accepted as isolated and pulse 8 would 
be rejected. This method provided an objective way to 
identify and quantify pulses that were influenced by neigh¬ 
bouring signals. 
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Fig. 2. The wavelet reconstruction of GRB 920513. (a) The GRB profile and wavelet fit with maxima (identified by 
the + symbol) and minima (identified by the x symbol). A total number of 33 pulses were identified above a threshold 
of 5cr. (b) The residuals between the wavelet fit and the actual GRB profile, (c)-(e) The decomposition of GRB 920513 
(BATSE trigger 1606) into its component wavelet scales. The wavelet transform for scales s = 2 2 - 4 (the finest scale, 
s = 2 1 , is not shown because it was dominated by noise), (f) The low frequency signal S^ 4 containing the remaining 
information for scales i > 4. 
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Fig. 3. A diagramatic representation of the algorithm for 
identifying pulses with appropriate significance. The error 
bars were calculated using the count rate in the particu¬ 
lar bin, multiplied by the threshold significance level, r a . 
The dashed circles highlight maxima and minima where 
the error bars overlap. Larger values of r a eliminate over¬ 
lapping regions that were more widely separated. With a 
slightly smaller choice of r CT in this example, pulses 6 and 
14 would have been accepted. 


2.4. Properties of the pulses 


The next task was to characterise the GRB profiles based 
on the properties of the constituent pulses. The follow¬ 
ing characteristics which had been studied previously 


(McBreen et al., 1994; Li & Fenimore, 1996; Hurley et al. 


199S| ) were investigated: the number of pulses per burst, 
N; the time intervals between pulses, AT; the pulse am¬ 
plitudes, C; the pulse area, A p \ the rise and fall times, 
t r and tf. and the pulse durations or full width at half 
maximum, FWHM; 

The total number of pulses in the sample of 319 GRBs 
was determined for a range of thresholds and varying 
isolation levels, tj (Fig. 4). The variation in isolation level 
has a much larger effect on the sample than the threshold, 
and caused a reduction in the number of pulses from over 
3000 at the isolation level of 20% to under 800 at the 
80% level. Fig. 4 also shows that the number of pulses 
falls quite quickly as r CT increases from 3 to 5. The initial 
rapid reduction in the number of pulses selected may be an 
indication of the removal of the small population of noise 
pulses remaining after the denoising process. The total 
number of pulses is not very sensitive to the threshold level 
in the region of 5er. The 319 GRBs are listed in Table 1, 
along with the total number of pulses above 5cr for each 
burst. Also included in the table are the number of isolated 
pulses at and above the 50% and 75% levels. 

In the analysis of pulse shapes, non-parametric meth¬ 
ods were used to estimate the various characteristics of 
the pulse profiles. This approach was chosen to make the 
conclusions more robust since no assumptions were made 
about the pulse shapes. Also since the measurements are 
made on the isolated pulses selected by the algorithm, the 
degree of isolation can be varied arbitrarily. If a particular 
measurement was sensitive to influence from surrounding 
pulses then the threshold tj was increased until the in¬ 


line 1 - 
3e+03 ■ 
2.5e+03 



Fig. 4. The number of pulses selected from the 319 GRBs 
as a function of the isolation level tj and the threshold 
t g . The plot also shows a projection of the data on to 
the tj , T(j plane and the contour levels are as given in the 
legend. 



No. of Pulses 

Fig. 5. The distribution of the number of pulses (N) per 
GRB. The shaded regions highlight the division into four 
categories, namely M, N, O and P with 1 < A < 2, 3 
< A < 12, 13 < A < 24 and N > 25 respectively. 


fluence was reduced sufficiently, with the proviso that the 
number of pulses in the sample remains statistically use¬ 
ful. 

The classification of pulses into isolated and non¬ 
isolated categories based on the algorithm allowed the 
measurement of characteristics of the temporal profile 
which are affected by neighbouring signals. The level of 
interference between pulses and the surrounding signal is 
dependent on the threshold at which the selection of these 
pulses is made. In fact a broad range of threshold levels 
were used to examine the time profiles and tj was typ¬ 
ically varied from 20% to 80%. It was decided based on 
these results to adopt pulses with t ct > 5er and tj > 50% 
for the main analysis of the pulse properties. 
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The pulse amplitude was measured as the maximum 
count in a 64 ms time interval after background subtrac¬ 
tion. The pulse area was measured using the sum of the 
background subtracted count rates starting at 5% of the 
height of the pulse above the left minimum to 5% of the 
pulse height above the right minimum on the falling edge 
of the profile. The starting point at 5% above the min¬ 
imum was chosen to eliminate contributions from back¬ 
ground noise for pulses with minima widely separated 
from the maximum. 

For similar reasons the rise time was measured from 
5% of the height of the pulse above the left minimum to 
95% of that height. The upper value of 95% ensures that 
the finishing point is robust against flat topped pulses and 
noise in the profile near the maximum. The fall time was 
measured in a similar way to the rise times i.e. from 95% 
to 5% of the pulse height above the right minimum. 

The duration of the individual pulses was measured us¬ 
ing the FWHM of the pulse. This approximation is valid 
only for well isolated pulses and tended to give poorer 
estimates for the true pulse width as the effect of neigh¬ 
bouring signals increased and the left and right minima of 
the pulse rose out of the background. 

2.5. Number of pulses per GRB 

The frequency distribution of the number of pulses (N) 
per GRB is given in Fig. 5. N has a range from 1 to 86 
with a peak at a value of 3, a median of 6 and only 10% 
of GRBs have N > 25. For convenience N is divided into 
four categories (Fig. 5). There are 34 GRBs in category P 
and only 7 have N > 50. Many of the timing studies on 
GRBs have concentrated on the categories with large N 
(e.g. Norris et al. 1996; Li & Fenimore 1996), which means 
that these analyses have focussed solely on the tail of the 
distribution shown in Fig. 5. 


3. Results 


3.1. The Lognormal distribution 

It has been shown previously that the lognormal distri¬ 
bution can adequately describe the properties of GRBs 
( iQuilligan et al., 1999 ). The lognormal distribution is gen¬ 
erated by statistical processes whose results depend on 
a prod uct of probabilities arisin g from a combination of 
events ( Aitehison fc Brown, 1957 ). A positive random vari¬ 
able X is said to be lognormally distributed if Y=log(X) 
is normally distributed with mean p and variance <r 2 . The 
probability density function is: 


/(x) = {^«P<-('»-M)V2Ox>0 (1) 

The median of the distribution occurs at x = e^. Many 
examples of lognormal distributions occur in nature, such 
as the propagation of a laser beam in a turbulent medium, 


the size of cumulus clouds in the atmosphere and ter¬ 
restrial lightning (Urnan, 1987; McBreen et al., 1994). In 
the case of terrestrial lightning, the durations, peak cur¬ 
rents, intervals between the strokes in the flashes, and the 
flash charges are all lognormally distributed. The statis¬ 
tical properties of strokes in lightning discharges in the 
Earths atmosphere have a close resemblance to the sta¬ 
tistical properties of pulses in GRBs. In both cases the 
result depends on a multiplicative process arising from a 
combination of independent events. 


3.2. Pulse analysis 
3.2.1. FWHM of pulses 

The distribution of the FWHM of the pulses from 319 
GRBs and the best lognormal fit are given in Fig 6 a. Fig. 
6 b shows the same data plotted as a cumulative percent 
on logarithmic probability paper such that a lognormal 
distribution yields a straight line. All pulses have r CT > 
5 and 17 > 50%. The distribution is very broad with a 
maximum at 0.7 s and half widths at 0.14 s and 3.5 s. 
The value of the reduced y 2 is 0.3 showing the data is 
well fit by the lognormal distribution. In obtaining x, 
the part of the fit below 0.1 s was not included because 
of the distortion caused by the limited time resolution 
that is apparent in Fig. 6 b. The value of the reduced x 2 
for the best lognormal fit as a function of T a and 77 are 
given in Figs. 6c and 6d. The fits are acceptable over most 
of the range with the largest departures occurring at the 
lowest values of and 77 because of the serious effects 
of pulse pile up. The lognormal distribution is a better 
fit when the effects of the overlapping pulses are reduced. 
The parameters of the best lognormal fit for r CT > 5 and 
ry > 50% are given in Table 2. 


3.2.2. Rise time and fall times of the pulses 

The distribution of the rise times of the pulses and the 
best lognormal fit are given in Fig. 7a. Fig. 7b shows the 
same data plotted as a cumulative percent and the large 
count in the first bin is due to the 64 ms resolution of 
BATSE. The parameters of the best lognormal fit to the 
broad distribution are listed in Table 2. The value of the 
reduced \ 2 = 1-3 show the data is compatible with the 
lognormal distribution. 

The distribution of the fall times of the pulses and 
the best lognormal fit are given in Fig. 8. The parameters 
of the best lognormal fit to the data are given in Table 
2. The fall times are also compatible with the lognormal 
distribution and have a wider range with a slightly larger 
mean than the rise times. 


3.2.3. Areas and amplitudes of the pulses 

The distributions of the amplitudes and areas of the pulses 
with T a > 5 and 77 > 50% from 250 GRBs, summed over 
two detectors, and the best lognormal fits are given in Fig. 
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Table 1. The BATSE trigger number of the GRBs used in the analysis. N is the total number of pulses in the GRB, 
N(50/75) are the numbers of pulses that are isolated at and above the 50% and 75% levels respectively. 


GRB 

N 

N ( 50 / 75 ) 

GRB 

N 

N ( 50 / 75 ) 

GRB 

N 

N ( 50 / 75 ) 

GRB 

N 

N ( 50 / 75 ) 

GRB 

N 

N ( 50 / 75 ) 

GRB 

N 

N ( 50 / 75 ) 

105 

4 

3/2 

1974 

2 

0/0 

2994 

36 

20/6 

3929 

3 

2/1 

6113 

5 

3/1 

7318 

12 

5/3 

109 

19 

10/2 

1997 

13 

8/2 

3001 

3 

3/2 

3930 

19 

6/1 

6124 

30 

14/5 

7329 

3 

0/0 

130 

11 

5/1 

2037 

6 

5/3 

3035 

21 

4/1 

3936 

5 

1/1 

6168 

2 

1/1 

7343 

11 

2/1 

143 

13 

6/3 

2053 

1 

1/1 

3039 

5 

3/1 

3954 

1 

1/1 

6198 

10 

1/1 

7360 

5 

4/1 

179 

4 

3/2 

2067 

8 

1/1 

3042 

19 

12/6 

4039 

33 

21/7 

6235 

6 

3/1 

7374 

1 

1/1 

219 

13 

4/3 

2080 

26 

12/2 

3057 

52 

3/1 

4048 

9 

3/3 

6242 

3 

3/3 

7429 

3 

1/1 

222 

6 

0/0 

2083 

2 

2/1 

3067 

8 

4/3 

4312 

4 

2/1 

6251 

3 

3/2 

7446 

1 

1/1 

249 

21 

4/1 

2090 

15 

7/2 

3105 

30 

19/9 

4368 

12 

2/2 

6266 

9 

2/1 

7464 

3 

2/1 

394 

25 

13/3 

2110 

13 

3/1 

3110 

13 

8/3 

4556 

3 

1/1 

6274 

18 

16/6 

7475 

2 

2/2 

451 

2 

2/2 

2138 

6 

2/2 

3115 

9 

6/3 

4701 

20 

8/4 

6303 

1 

1/1 

7477 

9 

3/2 

467 

3 

1/1 

2151 

3 

3/3 

3128 

36 

13/1 

4814 

2 

1/1 

6321 

3 

3/3 

7491 

40 

11/1 

469 

1 

1/1 

2156 

39 

8/3 

3138 

4 

2/1 

5080 

7 

7/3 

6329 

9 

1/1 

7503 

4 

1/1 

503 

1 

1/1 

2213 

5 

5/2 

3178 

6 

3/1 

5304 

14 

0/0 

6335 

2 

2/2 

7527 

2 

2/2 

543 

3 

2/1 

2228 

28 

12/5 

3227 

17 

5/3 

5389 

13 

7/3 

6336 

3 

0/0 

7529 

1 

1/1 

647 

3 

1/1 

2232 

14 

8/6 

3241 

15 

10/3 

5417 

2 

0/0 

6397 

1 

1/1 

7530 

2 

2/2 

660 

2 

2/1 

2316 

1 

1/1 

3245 

71 

24/6 

5419 

2 

1/1 

6400 

3 

3/3 

7549 

33 

6/2 

676 

4 

2/1 

2321 

11 

4/1 

3255 

12 

6/3 

5439 

4 

4/3 

6404 

13 

7/2 

7560 

12 

8/6 

678 

52 

21/6 

2329 

20 

3/1 

3269 

15 

13/7 

5447 

2 

1/1 

6413 

8 

7/3 

7575 

30 

10/4 

829 

3 

0/0 

2362 

2 

2/2 

3287 

10 

6/1 

5450 

13 

5/4 

6422 

2 

0/0 

7605 

9 

9/3 

841 

6 

6/3 

2367 

3 

3/2 

3290 

6 

4/3 

5451 

2 

0/0 

6436 

3 

3/3 

7607 

8 

8/3 

869 

10 

4/1 

2371 

2 

2/1 

3306 

16 

13/4 

5470 

7 

4/4 

6451 

3 

2/1 

7678 

37 

5/2 

907 

4 

3/2 

2387 

1 

1/1 

3330 

13 

5/3 

5473 

39 

27/19 

6453 

25 

5/1 

7688 

14 

9/3 

973 

3 

1/1 

2393 

2 

1/1 

3345 

4 

1/1 

5477 

12 

9/5 

6472 

47 

18/5 

7695 

27 

16/5 

999 

2 

2/2 

2431 

1 

1/1 

3351 

10 

7/5 

5486 

12 

1/1 

6525 

8 

6/5 

7711 

1 

1/1 

1025 

4 

2/1 

2436 

11 

8/2 

3408 

44 

17/5 

5489 

16 

8/4 

6528 

8 

3/1 

7766 

22 

15/6 

1085 

2 

0/0 

2446 

4 

3/2 

3415 

14 

8/4 

5512 

5 

1/1 

6560 

20 

16/9 

7770 

3 

2/1 

1122 

10 

3/1 

2450 

22 

7/4 

3436 

5 

4/1 

5523 

2 

2/2 

6576 

18 

11/5 

7775 

1 

1/1 

1141 

9 

0/0 

2533 

49 

9/1 

3458 

8 

2/2 

5526 

34 

34/14 

6587 

28 

9/4 

7781 

4 

2/2 

1157 

10 

8/2 

2537 

5 

2/1 

3480 

1 

1/1 

5530 

7 

5/1 

6593 

20 

13/2 

7788 

21 

14/6 

1159 

1 

1/1 

2586 

11 

7/7 

3481 

11 

1/1 

5548 

9 

8/4 

6621 

2 

2/2 

7845 

17 

15/11 

1190 

4 

2/2 

2611 

3 

2/2 

3488 

6 

6/5 

5563 

1 

1/1 

6629 

6 

4/3 

7858 

5 

0/0 

1204 

3 

3/3 

2628 

7 

5/2 

3489 

12 

2/2 

5567 

5 

1/1 

6630 

2 

1/1 

7884 

21 

15/6 

1385 

23 

4/1 

2700 

4 

3/1 

3491 

3 

1/1 

5568 

3 

0/0 

6665 

8 

2/1 

7906 

13 

7/5 

1419 

3 

1/1 

2736 

1 

1/1 

3512 

1 

1/1 

5572 

3 

3/3 

6672 

4 

3/1 

7929 

6 

4/1 

1425 

7 

5/1 

2790 

20 

10/1 

3516 

10 

7/2 

5575 

6 

2/1 

6683 

10 

10/3 

7954 

16 

14/13 

1440 

15 

9/2 

2793 

7 

4/3 

3523 

16 

0/0 

5591 

10 

7/4 

6694 

10 

6/1 

7969 

2 

1/1 

1443 

2 

2/2 

2797 

3 

1/1 

3569 

4 

2/2 

5593 

1 

1/1 

6764 

8 

3/1 

7987 

6 

2/1 

1468 

19 

17/6 

2798 

11 

1/1 

3593 

13 

6/1 

5601 

1 

1/1 

6814 

1 

1/1 

7994 

13 

1/1 

1533 

32 

20/12 

2799 

8 

6/6 

3598 

3 

3/2 

5614 

2 

2/2 

6816 

8 

7/3 

7998 

4 

3/2 

1541 

33 

H /1 

2812 

15 

12/9 

3634 

7 

7/3 

5621 

6 

3/1 

6824 

2 

2/2 

8008 

6 

3/3 

1578 

3 

0/0 

2831 

58 

12/2 

3648 

3 

3/1 

5628 

9 

4/1 

6904 

1 

1/1 

8019 

2 

2/1 

1606 

33 

15/4 

2852 

23 

5/1 

3649 

2 

2/2 

5644 

4 

1/1 

6930 

2 

2/2 

8022 

5 

5/4 

1609 

8 

2/1 

2855 

13 

2/1 

3658 

4 

2/1 

5654 

5 

3/2 

6963 

17 

10/2 

8030 

7 

6/3 

1625 

15 

2/1 

2856 

86 

19/3 

3663 

9 

7/3 

5704 

9 

8/4 

7028 

5 

3/2 

8050 

1 

1/1 

1652 

6 

1/1 

2889 

26 

12/4 

3765 

4 

2/2 

5711 

2 

2/1 

7113 

52 

8/1 

8063 

4 

4/3 

1663 

16 

0/0 

2891 

18 

15/4 

3776 

7 

1/1 

5725 

13 

12/6 

7170 

16 

6/2 

8087 

8 

7/3 

1664 

5 

1/1 

2894 

5 

1/1 

3788 

4 

3/1 

5726 

3 

3/1 

7172 

1 

1/1 

8098 

4 

2/2 

1676 

35 

21/7 

2913 

5 

2/2 

3860 

13 

9/3 

5773 

4 

1/1 

7185 

19 

16/7 

8099 

1 

1/1 

1683 

6 

3/1 

2919 

3 

3/2 

3866 

1 

1/1 

5867 

7 

6/3 

7240 

4 

4/4 

8111 

2 

2/1 

1709 

2 

0/0 

2929 

44 

16/5 

3870 

1 

1/1 

5955 

4 

2/2 

7247 

2 

1/1 




1711 

8 

2/1 

2953 

6 

2/2 

3891 

3 

2/1 

5989 

8 

7/5 

7255 

5 

1/1 




1815 

10 

5/2 

2958 

2 

2/2 

3893 

3 

1/1 

5995 

7 

1/1 

7290 

1 

1/1 




1883 

1 

1/1 

2984 

22 

16/6 

3905 

4 

4/1 

6090 

3 

1/1 

7301 

51 

13/7 




1886 

5 

0/0 

2988 

2 

2/2 

3912 

1 

1/1 

6100 

3 

0/0 

7305 

3 

3/3 





9 and Fig. 10. The distributions are very broad and the 
values of the best lognormal fits to the data are listed in 
Table 2. The lognormal distribution is compatible with 
the pulse areas and the amplitudes. 


3.2.4. T i me i nterva l s between pu l ses 


The distribution of the time intervals between the pulses 
is given in Fig. 11. The peak in the distribution occurred 
at about 1.0 s and was truncated at short time intervals 
by the 64 ms resolution of the data (Fig. 11a). A minimum 
time interval of 128 ms is required because two maxima 


must be separated by at least one time bin. There could 
be a large additional excess of pulses with separations be¬ 
low the 128 ms resolution of the data that were not re¬ 
solved. Other studies using different statistical methods 
from those employed here (Panaitescu et al., 1999; Lee 
et ah, 2000a||b|; |Spada et al., 2000D have also noted a deficit 


of time intervals below one second. A parent lognormal 
distribution of time intervals with parameters similar to 
the observed distribution was simulated and the time in¬ 
tervals between the pulses recorded with an accuracy of 
64 ms. The resulting distribution is given in Fig. 11a and 
the values of /i and cr for the parent distribution are given 


















F. Quilligan et al.: Temporal properties of GRBs 


9 



FWHM (sec) 



2.5 


3 
CT 
U1 

U 

0.5 


0 


1.5 


<V 

Jh 

cd 

3 1 

Cfi 

£ 

<J 0-5 


0 


Fig. 6. (a) The distributionpf the FWHM of pulses with 
t<t > 5 and rj > 50% and the best lognormal fit to the 
data, (b) The same data plotted as a cumulative percent 
such that a lognormal distribution yields a straight line. 
The large count in the first bin is due to the 64 ms time 
resolution, (c-d) The values of reduced y 2 for the best 
lognormal fit as a function of r G and 77 . 


(c) 


3 4 5 6 7 8 

Number of Sigma 


(d) 
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Isolation Level (%) 
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Fig. 7. (a) The distribution of the pulse rise times with 
t g > 5 and tj > 50% and the best lognormal fit to the 
data, (b) The same data plotted as a cumulative percent. 
The large count in the first bin is due to the 64 ms time 
resolution. 



t{ (sec) 


Fig. 8. The distribution of the fall times of the pulses with 
r,j > 5 and 17 > 50% and the best lognormal fit to the 
data. 

in Table 2. The measured distribution of time intervals in 
GRBs is consistent with the parent lognormal distribution 
provided a small (5%) excess of time intervals longer than 
15 seconds is not included. This excess is clearly visible in 
Fig. 11a. The time intervals greater than 15 s are plotted 
in Fig. lib. The data is well fit by a power law of slope 
- 1 . 2 . 
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Fig. 9. The distribution of the pulse amplitudes with r a > 
5 and tj > 50% and the best lognormal fit to the data. 



Pulse Areas (counts) 

Fig. 10. The distribution of the pulse areas with r CT > 5 
and tj > 50% and the best lognormal fit to the data. 

Table 2. The parameters of the best lognormal fit. The 
parameters are expressed as natural logarithms, with the 
width of the distributions shown in normal space. 


Property 


a 

x 2 

Width (±50%) 

FWHM 

-0.36 

1.37 

0.3 

0.14 - .3.5 

Rise time 

-0.44 

1.59 

1.3 

0.1 - 4.2 

Fall time 

-0.07 

1.59 

1.5 

0.14 - 6.1 

Pulse Amp. 

9.0 

1.12 

0.3 

2.2xl0 3 - 30x 10 3 

Area 

11.9 

1.2 

1.1 

35xl0 3 - 600 xlO 3 

Time Int.p] 

0.21 

1.03 

- 

0.37-4.14 

Peak Energyp| 

5.6 

0.58 

- 

137-535 


a See Sect. 

3.2.4 


6 See Sect. 

3.2.5 



3.2.5. Peak Energies 

The values of the pea k energy, E pea k, of a large sample 
of GRBs are given by Preece et al. (2000| ). There is an 
overlap of 77 GRBs with our sample, and the distribu¬ 
tion of the values of E pea k for each section of these bursts 
are given in Fig. 12. The distribution was well fit by a 
lognormal distribution with a small tail noticeable at low 




20 50 100 200 

Time Interval 

Fig. 11. (a) The measured distribution of time intervals 
between all the pulses and the best fit lognormal (dashed 
line) with allowance for the BATSE resolution of 64 ms. 
The excess of time intervals > 15 s is called the Pareto 
Levy tail, (b) The Pareto tail of the time intervals is well 
fit by a power law (dashed line) of slope ~ -1.2. 


values of E pea k. The values of /r and a are given in Table 2. 
The distribution of E pea k is noticeably narrower than that 
of the pulse parameters and spans a range of about 4 in 
width. 


3.3. Summary 


The major result of this part of the analysis is that the 
distributions of the rise times, fall times, FWHM, pulse 
amplitudes, pulse areas and the time intervals between 
the pulses are all very similar. The frequency distribu¬ 
tions are very broad and cover about three orders of mag¬ 
nitude and all are compatible with the lognormal distri¬ 
bution. Li fe Fenimore (1996| ) also showed that the pulse 
fluences and the time intervals between pulses are lognor- 
mally distributed for individual bursts in a small sample 
of bursts with more than 20 pulses. They also scaled the 
bursts to the same ^jl and a and showed that the sum¬ 
mation of all the peak fluences and time intervals looked 
lognormal, although no significance level was given for the 
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Fig. 12. The distribution of the peak energies and the best 
lognormal fit to the data (dashed line). 


result. No normalisation of pulse properties was applied to 
the GRBs in this analysis, because only the brightest 319 
bursts which had the best signal to noise ratio were anal¬ 
ysed. Fro m the data available for those bursts with known 
redshifts ( Atteia et al., 1999 ), there does not appear to 
be any dependable standard by which to scale the bursts 
because of the broad range of intrinsic luminosities and 
their comparatively small range of distances. Therefore, 
to avoid introducing further biases, and to use all of the 
pulse information available, the data were analysed with¬ 
out scaling. However, as a test of this process, the data 
were also scaled and the same analysis performed on the 
scaled data, and no significant differences were found be¬ 
tween the two data sets. 


3.4. Correlations between burst and pulse parameters 

It is important to determine how N relates to the other- 
parameters of the GRB. In Fig. 13 N is plotted versus 
burst duration (Tgo), total fluence and the median value of 
Epeak- Spearman rank order correlation co-efficients p and 
associated probabilities were obtained for the quantities in 
Fig. 13. The values are listed in Table 3 which also includes 
an additional range of burst parameters. The parameter 
C max is the maximum value of the peak amplitude in that 
burst. The high values of p show a strong correlation be¬ 
tween N and the total fluence, Tgo and E pea k. 

The values of p are not always uniformly distributed 
within each burst category. Tgo versus fluence is much 
better correlated for category N than either O or P. N 
versus C max and N versus HR are better correlated for 
category P than either N or O. 

The Spearman rank order correlation coefficients and 
probabilities were evaluated for isolated pulses with the 
range of pulse parameters given in Table 4. The pulse 
parameters are strongly correlated with each other. The 
pulse amplitude is negatively correlated with the pulse rise 
and fall times and FWHM. In general the correlations are 



Epeak ( keV ) 


Fig. 13. The number of pulses in a GRB as a function 
of (a) Tgo (b) total fluence and (c) the median value of 
the peak energy of the GRB. Quantisation occurs in the 
figures for low values of N. 


stronger for categories M and N than either O or P. The 
only significant exception to this trend is the pulse ampli¬ 
tude versus area which also has the highest values of p for 
categories O and P. 
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Table 3. Spearman rank order correlation coefficients be¬ 
tween a range of burst parameters 


Properties 

p 

Probability 

N VS. Tgo 

0.45 

5 x 10~ iy 

N vs. Total Fluence 

0.58 

3 x 10' 30 

N VS. C max 

0.27 

6 x 10- 7 

N vs. Hardness Ratio 

0.29 

1.2 x 10~ 7 

N VS. Epeak 

0.49 

7 x 10” 6 

Tgo vs. Total Fluence 

0.52 

6 x 10~ 24 

Tgo VS. C max 

-0.08 

0.16 

Tgo vs. Hardness Ratio 

0.11 

0.05 

Total Fluence vs. Maximum Pulse 

0.48 

1.3 x 10~ 15 

Total Fluence vs. Hardness Ratio (HR) 

0.56 

3 x 10~ 27 

C max VS. HR 

0.25 

4.4 x 10~ 6 


Table 4. Spearman rank order correlation coefficients be¬ 
tween a range of pulse parameters. 


Properties 

P 

Probability 

Rise Time vs. Fall Time 

0.64 

< 10 -48 

Rise Time vs. FWHM 

0.65 

< 10“ 48 

Rise Time vs. Pulse Area 

0.34 

2.9 x 10“ 34 

Rise Time vs. Pulse Amplitude 

-0.27 

1.9 x 10 -21 

Fall Time vs. FWHM 

0.70 

< 10' 48 

Fall Time vs Pulse Area 

0.42 

< 10~ 48 

Fall Time vs. Pulse Amplitude 

-0.22 

1.4 x 10 -14 

FWHM vs. Pulse Area 

0.44 

< 10“ 48 

FWHM vs. Pulse Amplitude 

-0.27 

2.7 x 10~ 21 

Pulse Area vs. Pulse Amplitude 

0.63 

< 10“ 48 

FWHM vs. Time Interval 

0.58 

00 

1 

O 

t-H 


3.5. Correlations between the time intervals between 
pulses and pulse amplitudes 

Spearman rank order correlation coefficients and probabil¬ 
ities were evaluated for the time intervals between pulses 
(AT) with T a > 5. The results are presented in Table 5 for 
two cases (1) the time intervals were not normalised and 
(2) the time intervals were normalised to Tgo- There is a 
good correlation between the time intervals in both cases 
that declines slowly with increase in the number of time 
intervals. The largest values of p occured in category N. 

The Spearman correlation coefficients were also eval¬ 
uated between pulse amplitudes and the results are given 
in Table 6 for two cases (1) the amplitudes were not nor¬ 
malised and (2) normalised to the largest amplitude pulse 
in the burst. The normalised pulse amplitudes are less 
strongly correlated over many pulses than the time inter¬ 
vals. These results were obtained for all pulses with r CT > 5 
and without satisfying any selection based on pulse isola¬ 
tion. The role of pulse pile-up has yet to be investigated. 

3.6. The properties of the pulses as a function of, N, 
the number of pulses in the GRB 

It was noticed early in th is analysis that pulse properties 
depended strongly on N ( |Quilligan et al., 2000| ). The me¬ 
dian value of the isolated pulse timing parameters were 


Table 5. Spearman rank order correlation coefficients p 
for time intervals between pulses. The two values for p 
and the probability are for unnormalised/normalised time 
intervals. 


Number of 
Intervals 

Total 

Number 

P 

Probability 

1 

2751 

0.42/0.56 

< 10~ 48 

2 

2499 

0.34/0.48 

< 10~ 48 

5 

1929 

0.24/0.37 

5 x 10- 26 / < 10 -48 

10 

1395 

0.20/0.29 

3 x 10 _13 /6 x 10 -27 

15 

890 

0.16/0.25 

3 x 10 _6 /4 x 10 -14 

20 

634 

0.10/0.23 

8 x 10 _3 /3 x 10" 9 

25 

459 

0.08/0.22 

7 x 10~ 2 /1 x 10 -6 

30 

322 

0.03/014 

3 x 10 2 /10" 2 


Table 6. Spearman rank order correlation coefficients 
p and associated probabilities for the pulse amplitudes. 
The two values for p and the probability are for unnor¬ 
malised /normalised amplitudes. 


Number of 
Amplitudes 

Total 

Number 

P 

Probability 

1 

3039 

0.72/0.57 

< 10“ 48 

3 

2499 

0.55/0.32 

< 10“ 48 

5 

2098 

0.52/0.24 

< 10 _48 /3 x 10“ 29 

7 

1777 

0.48/0.15 

< 10 —48 /6 x 10" 11 

9 

1510 

0.43/0.08 

< 10 _48 /2 x 10” 3 

10 

1395 

0.44/0.08 

< 10“ 48 /3 x 10~ 2 


determined for all GRBs with the same value of N. The 
median values of the rise time, fall time, FWHM and time 
interval between pulses are plotted versus N in Fig. 14(a- 
d). The largest value usually occurred for N = 1 or 2 and 
subsequently declined significantly as N increased. There 
are some values that are well removed from the general 
trend but they usually have a small number of pulses. 
The median values of the area and amplitude for isolated 
pulses are given in Fig. 15 (a-b). The trend is quite dif¬ 
ferent from Fig. 14. The amplitude is reasonably constant 
up to N ~ 35 with a clear increase for higher values of 
N. There is a similar but weaker trend for the pulse area 
which has the largest value at N = 1. 

The properties of the four categories of GRBs are sum¬ 
marised in Table 7. The median values of the pulse timing 
parameters all decrease by at least a factor of four from 
category M to P. In contrast the median values of Tgo, 
total fluence, hardness ratio and maximum pulse ampli¬ 
tude all increase significantly. The median variability, is 
defined as the number of pulses > 5 <r divided by the time 
the emission is > 5 cr, also increases from category M to 
P. 

3.7. First Half/Second Half analysis 

To study the evolution of the time profile as the GRB 
progresses, each GRB was divided into two, to include 
the pulses that occur before and after the strongest pulse 
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Table 7. The properties of the four categories of pulses in GRBs. The last five entries are for the 250 GRBs that were 
summed over two detectors. 


GRB Category 

M 

N 

O 

P 

Number of Pulses per GRB 

1-2 

3-12 

13-24 

25+ 

Number of GRBs 

67 

162 

56 

34 

Total number of pulses 

103 

981 

933 

1341 

Number of isolated pulses at 50% level 

83 

522 

476 

494 

Median T 90 (sec) 

18.1 

20.4 

45.7 

58.7 

Median Total Fluence (ergs/cm 2 ) 

8 . 8 xl 0 -6 

1.7xl0 -5 

4.2xl0 -5 

1 . 2 xl 0 -4 

Median hardness ratio (Chan |i|) 

3.4 

4.1 

6.5 

8.3 

Median Cmax (ph/cm 2 /sec) 

4.7 

6.1 

9.3 

12.6 

Median Variability 

0.09 

0.29 

0.39 

0.53 

Median Rise Time (sec) 

1.7 

0.8 

0.64 

0.45 

Median Fall Time (sec) 

5.2 

1.5 

1.0 

0.7 

Median FWHM (sec) 

1.8 

0.7 

0.64 

0.45 

Median Time Interval (sec) 

4.8 

1.9 

1.5 

1.0 

Number of GRBs 

55 

130 

38 

27 

Total number of pulses 

87 

778 

648 

1081 

Number of isolated pulses at 50% level 

70 

416 

319 

404 

Median Pulse Amplitude (Iso. pulses) (counts) 

13xl0 3 

7.3x10 s 

7.3 xlO 3 

9.5xl0 3 

Median Area (counts) 

560 xlO 3 

180 xlO 3 

140 xlO 3 

140 xlO 3 


Table 8. Summary of the first half/second half analysis of the three categories of the GRBs. 


GRB Category 

N 

O 

P 

Number of Pulses per GRB 

Total Number of Pulses (lst/2nd half) 

Total Number of Isolated Pulses (50%) (lst/2nd half) 
Total Number of Isolated Pulses (75%) (lst/2nd half) 

3-12 

404/415 

193/220 

79/126 

13-24 

384/493 

182/261 

77/101 

25+ 

679/628 

203/271 

71/83 

Median Rise Time (lst/2nd half) 

K-S Statistic/Probability 

1.02/0.70 

0.13/5% 

0.77/0.58 

0.12/9% 

0.51/0.45 

0 . 11 / 11 % 

Median Fall Time (lst/2nd half) 

K-S Statistic/Probability 

1.09/1.89 
0 . 20 /. 02 % 

0.90/1.22 

0.15/2% 

0.70/0.70 

0 . 11 / 10 % 

Median Asymmetry Ratio (lst/2nd half) 

K-S Statistic/Probability 

1 . 0 / 0 .4 
0.27/3 x 10" s % 

0.75/0.49 

.17/0.3% 

0.75/0.54 

0.14/1.5% 

Median FWHM (lst/2nd half) 

K-S Statistic/Probability 

0.77/0.67 

0.08/55% 

0.64/0.64 

0.08/52% 

0.51/0.45 

0 . 11 / 11 % 

Median Time Interval (lst/2nd half) 

K-S Statistic/Probability 

1.60/2.10 

0.13/2% 

1.41/1.54 

0.06/35% 

1 . 02 / 1.02 

0.05/30% 

Median Pulse Amplitude (x 10 3 )(lst/2nd half) 

K-S Statistic/Probability 

6 .2/5.4 
0.09/41% 

7.3/6 .8 
0.08/52% 

18.1/11.2 

0.09/25% 

Median Area (x 10 3 )(lst/2nd half) 

K-S Statistic/Probability 

151/127 

0.08/50% 

133/118 

0.07/59% 

152/130 

0.14/2% 


in the burst. Only GRBs with more than two pulses are 
included, resulting in a reduced total sample of 252. The 
first half (pre-main pulse) of the GRB was compared with 
the second half (post-main pulse). The bursts are also sub¬ 
divided into three categories. A summary of the properties 
of the GRBs used is given in Table 8. 

The first half/second half analysis was performed on 
the three timing parameters of the pulses, time intervals 
between the pulses, amplitude, area and the pulse asym¬ 
metry ratio which is defined as the ratio of the pulse rise 


time to the pulse fall time. The median values of the distri¬ 
butions in the three categories in the first half and second 
half analysis are given in Table 8 along with the results of 
the Kolmogorov-Smirnov (KS) tests. The KS probability 
is a measure of whether the two distributions (first/second 
half) are drawn from the same parent distribution. 

The first result is that the median values of the timing 
parameters of the pulses and the time intervals between 
the pulses all decrease by an average of 1.8, from the cat¬ 
egory N to P including the first half and second half of 
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Fig. 14. The median values of (a) rise time, (b) fall time, 
(c) FWHM and (d) the time intervals between pulses ver¬ 
sus the number of pulses 
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Fig. 15. The median values of (a) pulse area and (b) pulse 
amplitudes versus the number of pulses. 


the GRBs. In the case of the rise times, the median values 
of the distributions decrease from 1.02 to 0.51 in the first 
half and 0.70 to 0.45 in the second half. The trend in the 
median value of the pulse amplitude is in the opposite di¬ 
rection with larger pulses in category P than either O or 
N. 

In the first/second half analysis there is a trend in the 
three categories for the median rise time to be slower in 
the first half of the burst (1.02 versus 0.70 for category N). 
The difference could be caused by an additional clearing 
out effect at the start of the GRB. There is also a clear in¬ 
dication at the 0.02% level that the pulse fall time is faster 
in the first half than the second for category N (1.09 versus 
1.89 for category N) and this effect weakens for categories 
O and P. The median values of the pulse asymmetry ratio 
also show the most significant differences for category N 
where the median values are 1.0 and 0.4 for the first and 
second halves. The KS test gives good agreement between 
the first half and second half for the FWHM, time inter¬ 
vals between pulses, the pulse amplitudes and areas. The 
median values of the pulse amplitude and areas are how¬ 
ever larger in the first half than the second half for the 
three categories of GRBs. 


No. of Pulses 
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4. Discussion 

4.1. Pulse shapes in GRBs 

A complete study of the BATSE time profiles of the 
brightest 319 GRBs has been presented. The statistical 
analysis of the data reveal the ubiquitous nature of the log¬ 
normal distribution in GRB time profiles (Figs. 6-10). The 
means and variances of the best fit lognormal functions 
are given in Table 2. The results presented in Figs. 6-10 
are for the isolated pulses from the four GRB categories. 
The data in Table 7 and the first half/second analysis 
show that the median values of the pulse properties vary 
with N. Lognormal distributions also apply to the spectral 
properties of GRBs. The BATSE spectroscopic detectors 
revealed that the break energies in G RB spectra are com - 
patible with a lognormal distribution ( Preece et al., 2000|) . 
The FWHM of the distributions that describe the pulses 
are in the range 14-45 (Table 2) whereas the value for 
Epeak is only « 4 and extends from 137 keV to 535 keV. 
The unexpected narrowne ss of the E pea k d istribution is a 
major problem in GRBs (Brainerd, 2000). Furthermore, 
spectral fitting of 41 pulses in 26 GRBs showed that the 
spectral hardness parameter E pea k decays linearly with 
energy fluence and that the distribution of the decay con¬ 
stants is roughly lognormal ( |Crider et al., 199£ ). 

The Spearman rank order correlation coefficients were 
obtained for a range of burst parameters (Table 3). There 
is good agreement with the results of Lee et al. (2000a) 
who used a different method. N is strongly correlated with 
total fluence, T go , HR and E pea k (Fig. 13). N is an impor¬ 
tant quantity in determining GRB properties and provides 
the link between total fluence and duration since both in¬ 
crease with N. The correlation coefficients between the 
parameters that describe the pulses are given in Table 4. 
The rise and fall times, FWHM and area of the pulses are 
highly correlated. There is a high probability that a pulse 
with a fast rise time will also have a fast fall time and 
a short FWHM. The pulse amplitude is negatively corre¬ 
lated with the rise and fall times and FWHM. The anticor¬ 
relation between the pulse amplitude and pulse width has 


been observed in other studies (Lee et al., 2000a; Fenimore 


fc Ramirez-Ruiz, 2001). The pulse width is also a func- 
tion of energy and varies as E~°- 45 (Fenimore et al., 1995) 


and th is effect has been attributed to synchrotron radi¬ 
ation ( Piran, 1999 ). The FWHM is strongly correlated 
with the preceding and subsequent time interval between 
pulses (Table 4). This result is in agreement with Nakar 
& Piran (2001) and the prediction of the internal shock 
model. A further comparison of the pulse timing param¬ 
eters, the energy dependence and the spectral lag may 
reveal further interesting constraints on the emission pro¬ 
cess (Norris et al., 2001; Daigne & Mochkovitch, 2001). 
The unique properties of the pulses in GRBs have been 
summarised by McBreen et al.,(2002). 

The lognormal distribution arises from a statistical 
process whose result depends on a product of probabil¬ 
ities arising from a combination of independent events. It 


therefore identifies the statistical process but not the com¬ 
bination of events that lead to the formation of the pulse 
shape and the peak energy. In the internal shock model 
the main factors contributing to the pulse shape include 
(Rees & Meszaros, 1994; Piran, 1991). 


1. the Lorentz factors, masses and thickness of the inter¬ 
acting shells 

2. the distance from the central engine and curvature 
where the collisions occur 

3. the energy conversion in the shock, the nature of the 
magnetic field and particle acceleration process 

4. the synchrotron radiation mechanism possibly modi¬ 
fied by self-absorption and pair production 

5. the time scale for energy loss by the particles. 

The resulting properties of the pulses in GRBs depend 
on a combination of many events and hence it is not sur¬ 
prising that the lognormal distribution gives an elegant 
description of all their properties. Random multiplicative 
processes abound in a variety of natural phenomena and a 
good example is the statistical description of the strokes in 
flashes of lightning. Almost all the properties of the strokes 
in the flashes and the flashes themselves are well described 


by lo gnormal distributions ( Uman, 1987 ; McBreen et al. 

1991 ). 


4.2. Time intervals between pulses in GRBs 

It is often found that distributions that seem to be lognor¬ 
mal over a wide range change to an inverse power law dis¬ 
tribution for the last few percent. An amplification model 
has been used to characterise the transition from a log¬ 
normal distributi on to a power law that is ofte n called a 
Pareto-Levy tail ( Montroll fe Shlesinger, 1982 ). The dis¬ 
tribution of time intervals conform to the lognormal distri¬ 
bution over most of the range with the exception of about 
5% of the time intervals longer than about 15 s (Fig. 11). 
The Pareto-Levy tail of time intervals have an amplifica¬ 
tion process that is not available to most time intervals. 

The origin of the nonrandom distribution of time in¬ 
tervals between pulses is an important clue to the GRB 
process. In the internal shock model there is almost a one 
to one correspondence between the emission of shells and 


pulses resulting from the collisions of shells (Kobayashi 
|et al., 19~97). Hence the time intervals between pulses is 


an almost direct measure of the activity of the central en¬ 
gine. The temporal behaviour of soft gamma-ray repeaters 
and young pulsars provide additional context in which to 
view the results of GRB time profiles. The time inter¬ 
vals between about 30 microglitches in the Vela Pulsar are 
consistent with a lognormal distribution with a mean of 
50 days ( |Hurley et al., 1994 , 1995; Cordes et al., 1988). 
The amount of energy involved in the microglitches is 
about 10 38 ergs. The macroglitches in the Vela Pulsar are 
about a thousand times more powerful but occur too in¬ 
frequently to determine the distribution of time intervals 
but they have a wide range and do not seem inconsis¬ 
tent with the lognormal distribution. More energetic out- 
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bursts have been recorded from SGR sources. The two 
most energetic events released about 5 x 10 44 ergs in 7 -rays 


from SGR 0526-66 on 5 March 1979 (Mazets et al., 1979) 
and about 10 43 ergs from SGR 1900+14 on 27 August 
1998 ( Hurley et al., 1999| ). The SGR sources also gener¬ 
ate a large number of smaller outbursts, and it has been 
shown that the time intervals between outbursts are dis¬ 


GRB models leave open many possibilities to account 
for the Pareto-Levy tail of long time intervals. The excess 
of long time intervals have been noted in other studies 
(Ramirez-Ruiz & Merloni, 2001; Nakar & Piran, 2001). 
The properties of the GRBs with long time intervals will 
be covered in a separate publication. 


tributed lognormally (Hurley et al., 1994; Gogiis et al. 


2000[)~|Hence lognormally distributed time intervals be- 4.3. Numbers of pulses and jets in GRBs 


tween outbursts and glitches are characteristic features of 
SGR sources and neutron star microglitches. It is widely 
accepted that these sources are rotating neutron stars with 
high magnetic fields. It is not unreasonable to argue that 
the co upled effects of rapid rotation and intense magnetic 
fields ( Kluzniak fe Ruderman, 1998 ) are also involved in 
powering GRBs since the time intervals between pulses 
are also consistent with a lognormal distribution. 

The possibility that a rapidly rotating neutron star 
with a surface magnetic field o f ~ 10 15 Ga uss could power 
a GRB has been suggested (Usov, 1992). Once formed 


such a neutron star could lose its rotational energy catas¬ 
trophically on a time scale of seconds. The rotation of the 
star decelerates because of the applied torques. Powerful 
transient fields may also occur in the merger of two neu¬ 
tron stars or a neutron star and a black hole. The en¬ 
ergy stored in differential rotation of the collapsed object 
would be released in sub-bursts as toroidal magnetic fields 
are repeatedly wound up to ~ 10 17 Gauss (Kluzniak & 
Ruderman, 1998). The emergence of a toroid is accompa¬ 


nied by huge spin down torques, the reconnection of new 
surface magnetic fields and rapid release of a sub-burst 
of energy of about 10 51 ergs. The release of rotational 
energy in repeated sub-bursts could power the GRB. If 
the differentially rotating compact object forms a torus 
about a spinning black hole either in a merger or core 
collapse of a massive star, energy can be extracted by 
the magnetic field th at threads the torus and the black 
hole ( Meszaros, 2001 ). As the torus builds up and ejects 
its magnetic toroids, the differential rotation of the torus 
could be maintained by the spin of the black hole.The 
models of GRBs with the coupled effects of rapid rotation 
and ultra intense magnetic fields are particularly attrac¬ 
tive because the time intervals between pulses in GRBs are 
distributed lognormally and follow the pattern observed in 
non-catastrophic events in SGRs and pulsars. 

The time intervals between the pulses are correlated 
with each other and the correlation decreases slowly with 
increase in the number of time intervals (Table 5). This 
effect had previously been observed in a small sample of 
GRBs ( Nakar fe Piran, 200l| ) and attributed to the inter¬ 
nal shock model. In addition the pulse amplitudes are also 
correlated with each other and this effect decreases more 
rapidly than the time intervals between pulses (Table 6 ). 
Similar correlations have been found between the pulse 
amplit udes and also time int ervals between pulses in short 
GRBs ( |McBreen et al., 2001 ). In the internal shock model, 
these correlations originate in the central engine and pro¬ 
vide strong constraints on any viable model of GRBs. 


A detailed comparison has been made between the dis¬ 
tributions of the properties of the pulses in the first half 
and second half for three categories of GRBs. There are 
no statistically significant differences between the median 
values of the time intervals between pulses, pulse ampli¬ 
tude, areas and FWHM in the first half and second half 
of GRBs (Table 8 ). There are two trends in the pulse rise 
times and fall times that should be noted: 1 ) the median 
rise time is slower in the first half for the three categories 
of bursts (Table 8 ) and 2) the median fall time is faster in 
the first half for categories N and O. The combination of 
slower rise times and faster fall times gives a pulse asym¬ 
metry ratio with a significant difference between the first 
half and second half for category N and at a reduced sig¬ 
nificance level for O and P (Table 8 ). The effect could be 
caused by a clearing out process such as additional baryon 
loading or Compton drag in the first half of the bursts 
with small number of pulses. These results are also com¬ 
patible with the constancy of the pulse widths observed 
by Ramirez-R.uiz fe Fenimore (2000| ) using a peak aligned 
profile method on a small sample of GRBs with more than 
20 pulses. In the internal shock model, the rapid variabil¬ 
ity in GRB time profiles is due to emission from multi¬ 


ple shocks in a relativistic wind (Piran, 1999; Panaitescu 
|et al., 1999 ; Downes et al., 2001). The temporal position 
of the pulse is unconnected to the collision parameters and 
in this way the litt le or no evolution of the pulses in GRBs 
can be explained ( Fenimore fe Ramirez-Ruiz, 2001 ). The 
rise time and fall time may be determined by the hydro¬ 
static time ssd/c and the angular spreading time ~D/c, 
where d and D are the w idth and separation of the shells 
( Kobayashi et al., 2001 ). The observed evolution of the 
pulses requires the shells to be narrower and farther apart 
later in the GRB. This prediction is in agreement with the 
data because the time intervals between pulses are longer 
in the second half of the burst (Table 8 ). 

However, it is evident from Figs. 14 and Tables 7 and 8 
that as the number of pulses in a burst increases, the me¬ 
dian values of the rise and fall times, FWHM and time 
intervals all decrease and by about the same amount. 
The GRBs with more pulses also have on average signifi¬ 
cantly longer durations, higher fluences and hardness ra¬ 
tios (Table 7). The variability index of a GRB was taken 
to be the number of pulses > 5cr divided by the time 
the GRB emission was also > 5tr. The median values are 
given in Table 9 for the four GRB categories. The GRBs 
with more pulses have a higher variability index. These re¬ 
sults provide an interesting interpretation of the two cor- 
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relations that have been reported for GRBs with known 
redsliift: (1) the more luminous GRBs to be more vari¬ 
able (Fenimore & Ramirez-Ruiz, 2001) and (2) there is 
an anticorrelation between the arrival times of high en¬ 
ergy and low energy pulses in GRBs (Norris et al., 2000; 


Salmonson, 2000). Recently Schaefer et al. (2001) showed 


that there is a relationship between the variability and 
spectral lag (loka & Nakamura, 2001). There is a good 
correlation betw een the values of the variability o btained 
here and those of Fenimore fe Ramirez-Ruiz (2001 ). GRBs 
with higher values of HR have lower values of < V/ V max > 
i mplying they a re a more distant and luminous population 
( [Schmidt, 2001 ). 

Our knowledge on the shape of the emitting region in 
GRBs is restricted because, due to relativistic beaming, 
only a small portion of angular size ~ F _1 is visible to 
the observer. Thus the observer is unable to distinguish a 
sphere from a jet as long as T > 9~ l where 9 is the radius 
of the opening angle of jet ( |Rhoads, 1997 ). However as 
the source continues its rapid expansion, T will decrease, 
and when F < fG 1 there will be a marked decrease in 
the observed flux. The steep time dependence of the af¬ 
terglow emission, sometimes with changes in the slope of 
the spectrum, and radio emission have been widely inter¬ 


preted as evidence for emission from jets (Castro-Tirado, 


2001; Meszaros, 2001; Frail et al., 2001). The GRBs with 


more pulses appear to have higher values of the Lorentz 
factor T. The higher values of T may come from a more 
efficient and more active central engine. In one variation of 
the internal shock model, it was assumed that the degree 
of collimation of the jet depended on the mass M at the 
explosion (Kobayashi et al., 2001). A wide jet involves a 
large mass that results in a flow with a lower F. The pulse 
properties depend strongly on r m i n ,r max and the radius 
of the photosphere R±. GRBs with faster pulses originate 
in collisions above R± whereas GRBs with slower pulses 
have smaller values of T and some collisions below the 
photosphere. While this homogeneous model may explain 
pulse properties in GRBs, the strong possibility of inho¬ 
mogeneous jets with a variable F should also be examined. 

Baryon loading can be a major problem in GRB mod- 
els and severely limit the attainable value of T ( Rees, 1999 ; 
Meszaros et al., 1998| ; Salmonson, 2000 ). There maybe a 
broad range of T’s in the outflow with the highest value 
occurring close to the rotation axis where the baryon con¬ 
tamination should be at a minimum. At larger angles from 
the axis, there may be an increasing degree of contamina¬ 
tion with a corresponding drop in F. The outcome of a 
collapse in a massive star whose iron core collapsed to a 


black hole have been computed (MacFadyen & Woosley 
1999[)7]rhe resulting jet that drives out through the star is 
probably powered by a MHD process which can in princi¬ 
ple convert a large portion of the binding energy at the last 
stable orbit into jet energy. The large amount of energy 
dumped into the natural funnel-shaped channel creates a 
highly collimated jet, focused into a small region of the 
sky. The largest value of F occurs on axis and decreases 
with increasing 0 because the material coming at the ob¬ 


server has less energy at larger angles. The emission is 
still beamed into an angle T _1 but in this inhomogeneous 
model the angle varies across the opening angle of the jet 
(Rossi,Lazzati & Rees, 2001). In this situation the prop¬ 
erties of the pulses in GRBs can be influenced by the jet. 
The BATSE sample of the brightest GRBs should contain 
a range of angles within the jet and hence different values 
of r. In this context it is reasonable to identify the com¬ 
plex GRBs with more pulses and higher values of T with 
angles near the axis of jet. The GRBs that are viewed at 
larger angles from the jet axis have on average lower val¬ 
ues of r, and develop at the greater distances from the 
central engine and should have slower pulses. 

In this context it is interesting to note that the pulse 
evolution consisting of slower rise times and faster fall 
times in the first half, is more pronounced for GRBs in 
category N than either O or P. In a jet model with a 
variable F, the GRBs in category N would be on average 
farther from the axis and more sensitive to a clearing out 
effect such a s additional bary on loading or Compton drag 
( Rees, 1999 ; Ghisellini, 200l| ) in the initial phase of the 
GRB. 

The steep and variable slope of the decay of GRB af¬ 
terglows have been widely interpreted as evidence for jets 
in GRBs (Meszaros, 2001; Castro-Tirado et al., 2001). If 
the axis of the jet is pointed close to the observer, the GRB 
will be intense and the afterglow should contain evidence 
for good alignment. It is interesting that the two bright¬ 
est GRBs detected by WFC on BeppoSAX also were the 
best aligned. The recent detection of a bright GRB with 
a fluence of 10 -4 ergs cm -2 also had a very bright af¬ 
terglow (Castro-Tirado et al., 2001). However many more 
GRBs and afterglows are required to verify the existence 
of a pa ttern between the strongest GRBs and their after¬ 
glows ( Frail et al., 2001 ). The distribution of the number 
of pulses per GRB (Fig. 5) may broadly represent the 
beaming by the jet because bursts with large numbers of 
pulses and higher variability (Table 7) may be close to the 
axis and bursts with smaller numbers of pulses and less 
variability further off-axis. 


5. Conclusions 

The properties of the brightest 319 GRBs in the BATSE 
current catalogue have been analysed.The automatic pulse 
selection process detected more than 3300 pulses. The 
distributions of pulse rise and fall times, FWHM, areas, 
amplitudes and time intervals between pulses are rea¬ 
sonably consistent with the lognormal distribution. GRB 
pulse profiles can be elegantly described by a small num¬ 
ber of parameters that may be very useful for simula¬ 
tions. The lognormal distribution depends on the product 
of probabilities arising from a combination of indepen¬ 
dent events and these conditions must therefore apply to 
the generation of the temporal and spectral properties of 
GRB pulses. A wide range of burst parameters and also 
pulse parameters were correlated and the results follow the 
trend expected from the internal shock model. The pulse 
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amplitude is strongly anticorrelated with the other pulse 
timing parameters. The time intervals between pulses and 
pulse amplitudes are correlated with each other. 

A comprehensive analysis has been performed between 
the first half and second half of GRBs in three categories 
defined in terms of N. No major differences were found 
between the distribution of pulse properties between the 
first half and second half of the GRBs. There is a strong 
tendency for pulses to have slower rise times and faster fall 
times in the first half of the burst. This trend is stronger in 
GRBs with small numbers of pulses. The pulse timing pa¬ 
rameters and time intervals all decrease with increase in N. 
These results seem to be compatible with jet models with 
either a T that varies with the opening angle or is constant 
and varies with the mass. If T varies with the opening an¬ 
gle of the jet, the GRBs with higher values of T and greater 
variability are observed close to the axis of the jet while 
GRBs with smaller number of pulses and less variability 
are observed at larger angles from the jet. Jets with val¬ 
ues of T that vary with angle or with mass may explain 
the luminosity-variability correlation and the luminosity- 
energy lag correlation in GRBs with known redshift. 

This study of the number of pulses in GRBs and their 
time structure provides strong evidence for rotation pow¬ 
ered systems with intense magnetic fields and the added 
complexity of a jet. These results can be well interpreted 
by internal shocks in the framework of theoretical models 
for the formation of black holes and subsequent jet forma¬ 
tion. 
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